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Dielectric  sensing  in  e-near-zero  narrow  waveguide  channels 

Andrea  Alu  and  Nader  Engheta* 

Department  of  Electrical  and  Systems  Engineering,  University  of  Pennsylvania,  Philadelphia,  Pennsylvania  19104,  USA 

(Received  28  May  2008;  published  3  July  2008) 

We  exploit  here  the  dramatic  field  enhancement  caused  by  energy  squeezing  and  tunneling  (i.e.,  “supercou¬ 
pling”)  in  metamaterial-inspired  ultranarrow  waveguide  channels  with  near-zero  effective  permittivity  in  order 
to  sense  small  permittivity  variations  in  a  tiny  object.  The  supercoupling  effect  is  accurately  modeled  analyti¬ 
cally  and  closed-form  expressions  are  derived  to  describe  the  presence  of  defects  or  permittivity  perturbations 
along  the  channel.  Applications  for  tailoring  its  pass-band  frequency  and  for  high-<2  sensing  are  proposed  at 
microwave  frequencies. 

DOI:  10.11 03/Phy sRe vB .78 .045 102  PACS  number(s):  41.20.Jb,  42.82.Et,  52.40.Db,  52.40.Fd 


I.  INTRODUCTION 

The  possibility  of  squeezing  energy  through  subwave¬ 
length  apertures  and  channels  has  been  recently  investigated 
in  several  systems  and  setups  involving  plasmonic  materials 
and  metamaterials  (see,  e.g..  Refs.  1-3).  Somewhat  analo¬ 
gous  to  these  findings,  the  idea  of  metamaterial  supercou¬ 
pling,  i.e.,  energy  squeezing  and  tunneling  through  subwave¬ 
length  narrow  waveguide  channels  of  arbitrary  length  and 
narrow  height  filled  with  e-near-zero  (ENZ)  metamaterials, 
has  been  put  forward  in  a  recent  letter.4  Further  theoretical 
results  and  experimental  realization  of  this  anomalous  phe¬ 
nomenon  have  been  obtained, 5-7  demonstrating  how  it  is 
possible  to  induce  such  tunneling  in  various  waveguide  con¬ 
figurations  at  microwave  frequencies.  In  particular,  we  have 
reported  in  Refs.  5  and  6  how  the  same  tunneling  phenom¬ 
enon  may  be  achieved  inside  a  simple  hollow  rectangular 
channel  of  narrow  height  that  operates  near  the  cutoff  fre¬ 
quency  of  its  dominant  TE10  mode,  thus  making  use  of  its 
natural  dispersion  to  mimic  the  response  of  ENZ  materials.8-9 
This  may  lead  to  resonant  transmission  and  vanishingly 
small  phase  delay  through  arbitrarily  narrow  waveguide 
channels,  independent  of  their  total  length  and  geometry. 

As  reported  and  discussed  in  Refs.  4-7,  this  energy 
squeezing  is  associated  with  a  huge  increase  in  the  magni¬ 
tude  of  the  electric  field  inside  the  ENZ  channels.  Intuitively, 
one  can  expect  that  such  high-intensity  field  may  ensure  high 
sensitivity  to  small  variations  in  material  parameters,  which 
may  then  be  exploited  for  material  sensing  applications.  Here 
we  analyze  in  detail  the  perturbative  effects  of  a  small  defect 
(e.g.,  a  small  change  in  material  permittivity)  or  of  a  small 
cavity  carved  inside  the  metamaterial  channel,  with  particu¬ 
lar  attention  to  the  resulting  changes  in  the  channel  tunneling 
properties.  In  particular,  in  the  case  of  high-quality-factor 
(i.e.,  high-Q)  tunneling,  it  is  shown  how  the  high-intensity 
electric  field  in  the  channel  may  be  successfully  used  for 
accurate  sensing  applications  and  for  fine  tuning  of  the  tun¬ 
neling  mechanism. 

II.  THEORETICAL  ANALYSIS  AND  NUMERICAL 
RESULTS 

Consider  the  geometry  of  Fig.  1,  consisting  of  an  ultrana¬ 
rrow  rectangular  channel,  of  height  c/ch  and  length  L,  con¬ 


necting  the  two  sections  of  a  waveguide  of  height  a  S>  ach.  All 
the  waveguide  sections  have  the  same  width  b>a,  which 
tailors  the  dispersion  properties  of  the  fundamental  TE10 
mode.  In  particular,  the  TE10  propagation  constant  in  each 
section  of  the  waveguide  may  be  regarded  as  that  of  a  TEM 
wave  propagating  in  an  unbounded  medium  or  in  a  parallel- 
plate  waveguide  with  effective  constitutive  parameters,8 

<W«b=  e-c2/(4f2b2), 

/A  ff=/An  (1) 

where  e  is  the  permittivity  of  the  homogeneous  material  fill¬ 
ing  the  waveguide  section,  c  is  the  velocity  of  light  in  free 
space,  and  an  e~'2"7>  time  convention  has  been  assumed.  At 
the  cutoff  frequency  of  this  mode,  i.e.,  when  b=kl 2  with  A. 
being  the  wavelength  in  the  filling  material  with  permittivity 
e,  the  effective  permittivity  (1)  “experienced”  by  the  mode  is 
zero,  consistent  with  the  infinite  phase  velocity  of  a  cutoff 
mode  or  of  a  TEM  wave  in  an  ENZ  metamaterial.  This  direct 
analogy  may  be  exploited  to  realize  metamaterial  effects  in 
several  geometries  (see,  e.g.,  Refs.  8  and  9).  Metamaterials 
realized  with  this  technique  are  inherently  robust  to  losses, 
since  they  avoid  the  use  of  small  inclusions  typical  of  con- 


FIG.  1.  (Color  online)  Geometry  of  the  structure  proposed  as  a 
permittivity  sensor.  The  two  input  and  output  rectangular 
waveguides  are  connected  via  an  ultranarrow  channel  (shown  in 
blue).  In  the  middle,  an  object  (i.e.,  a  cavity  region)  with  permittiv¬ 
ity  ecav  different  from  that  of  the  channel  is  inserted.  The  structure  is 
all  surrounded  by  conducting  walls. 
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ventional  metamaterial  technology,  which  usually  lead  to  ab¬ 
sorption  and  radiation  losses  associated  with  technological 
imperfections  and  disorder. 

As  we  have  done  in  Ref.  6,  this  technique  may  be  applied 
to  experimentally  verify  the  anomalous  properties  of  effec¬ 
tive  ENZ  channels  for  resonant  tunneling  and  zero-phase  en¬ 
hanced  transmission,  in  principle  independent  of  their  geom¬ 
etry  and  length.  In  the  geometry  of  Fig.  1  we  have  also 
considered  two  narrow  sections  of  length  ach  attached  to  the 
channel,  with  the  aim  of  “matching”  the  incident  mode  into 
the  narrow  channel,  forming  a  characteristic  U  shape  for  the 
ENZ  region  (blue  region  in  Fig.  1). 

In  this  way,  we  may  obtain  the  ENZ-related  supercoupling 
effects4-6  using  conventional  materials,  e.g.,  a  dielectric  or 
even  free  space  with  ech=  e0,  filling  the  U  channel.  Since  the 
ENZ-based  tunneling  is  expected  to  occur  at  the  cutoff  fre¬ 
quency  of  the  narrow  channel,  it  is  sufficient  to  have  the 
remaining  part  of  the  waveguide  filled  with  a  regular  dielec¬ 
tric  with  e>  ech  to  ensure  propagation  in  these  outer  sections 
of  the  waveguide  at  this  frequency.  The  phase  delay  between 
the  entrance  and  the  exit  of  the  narrow  channel  is  expected  to 
remain  very  small,  independent  of  the  channel  length,  due  to 
the  infinite  phase  velocity  of  the  mode  near  its  cutoff,  result¬ 
ing  in  a  uniform  and  strongly  enhanced  electric  field  along 
the  channel. 

Now  we  consider  a  small  cavity  region  with  length  Lcav 
and  permittivity  6cav  A  ech  in  the  channel  (Fig.  1).  The  pres¬ 
ence  of  this  cavity  region  with  a  permittivity  different  from 
the  rest  of  the  channel  clearly  perturbs  the  tunneling  effect. 
In  Ref.  6,  we  have  studied  the  effects  of  an  air  cavity  in  the 
ENZ  narrow  channel  on  the  electric-  and  magnetic-field  dis¬ 
tributions  in  the  channel.  Here  we  concentrate  on  how  and  to 
what  extent  the  cavity  with  ecav  A  ech  may  alter  and  tune  the 
tunneling  and  supercoupling  aspects  of  the  ENZ  channel.  As 
it  will  be  shown  in  the  following,  this  analysis  may  provide 
us  with  quantitative  measures  of  the  perturbation  caused  by 
the  presence  of  a  small  object  (here  a  cavity)  of  permittivity 
6cav.  This  may  directly  lead  us  to  useful  applications  of  this 
setup  for  sensing  dielectric  properties  of  an  object  under  test. 

The  overall  system  of  Fig.  1  may  be  analyzed  using  a 
transmission-line  (TL)  model,  consistent  with  the  one  pre¬ 
sented  in  Ref.  7  for  a  different  geometry  of  the  supercoupler, 
since  each  waveguide  section  supports  the  same  dominant 
TE10  mode.10  In  particular,  each  section  is  described  by  its 
characteristic  impedance  r/=  'J/XqI  ecfr  and  wave  number  [3 
=  2t r/V/U.0eefT,  and  the  steps  in  the  U  channel  cross  section 
may  be  modeled  with  voltage  transformers  with  a  transform 
ratio  s  =  achla  and  parallel  reactive  loads  that  take  into  ac¬ 
count  the  localized  higher-order  evanescent  mode  excitation. 
In  this  specific  geometry,  since  the  steps  are  embedded  in  an 
effective  ENZ  region  (the  U  channel)  at  the  resonant  fre¬ 
quency,  these  parallel  loads  may  be  neglected,  since  the  con¬ 


dition  eeff—  0  implies  an  extremely  high  impedance  for  these 
loads.11  Similarly,  due  to  the  condition  ach<827r//3ch,  it  is 
possible  to  neglect  the  presence  of  the  small  lateral  transition 
regions  of  the  U  in  the  TL  model.  We  note  that  this  same 
model  (and  the  following  discussions)  applies  for  a  parallel- 
plate  waveguide  when  filled  with  metamaterials  with  permit¬ 
tivity  eeff,  consistent  with  the  geometry  in  Ref.  4. 

In  the  absence  of  the  cavity,  or  when  ecav=  ech,  the  input 
impedance  at  the  entrance  of  the  channel  may  be  written  as 

2s2vlhVout- isych(sych-  ?7out)(A77ch+  ’7out)sin(2/?chL) 

s2  rflh  +  7] 2ut  +  (s  7fch  -  r/0J ( s  r/ch  +  ?7out)cos(2/3chL) 

(2) 

where  the  subscripts  ch  and  out  refer  to  the  channel  and  the 
outside  region,  respectively.  Total  transmission  is  obtained 
when  Zin=  ?7out,  i.e.,  either  when  L=2ttI  (3ch  or  when 


s 

The  first  solution  corresponds  to  classic  Fabry-Perot  reso¬ 
nances  of  the  channel,  very  sensitive  to  the  length  of  the 
channel  and  the  frequency  of  operation,  whereas  condition 
(3)  is  a  unique  matching  condition  obtainable  only  near  the 
channel  cutoff  frequency,  which  is  interestingly  independent 
of  its  geometry  (and  in  particular  its  length).  The  strong  mis¬ 
match  at  the  entrance  and  exit  faces  of  the  channel,  due  to 
the  difference  in  heights,  is  compensated  by  the  huge  in¬ 
crease  in  its  characteristic  impedance  for  frequencies  right 
above  the  cutoff.  In  the  limit  of  ach— >0,  i.e.,  ,v  — >  0 ,  condition 
(3)  requires  the  system  to  operate  at  the  frequency  for  which 
6eri  =  0,  i.e.,  at  the  channel’s  cutoff  frequency,  independent  of 
the  channel  length.  This  is  consistent  with  the  theoretical 
results  obtained  in  this  limiting  case  in  Ref.  4  using  a  differ¬ 
ent  approach.  For  finite  but  narrow  enough  channels  (s  — 0), 
Eq.  (3)  suggests  operation  at  frequencies  for  which  eeff— 0+. 
It  should  be  mentioned  that  the  steep  variation  of  7)ch  right 
above  the  cutoff  frequency  implies  that  the  tunneling  fre¬ 
quency  indeed  occurs  very  close  to  the  cutoff  frequency,  ir¬ 
respective  of  variations  in  ach  (while  ach  <  a),  consistent  with 
the  theoretical  and  experimental  results  in  Refs.  4-7.  The 
presence  of  the  parallel  loads,  representing  the  E-plane  dis¬ 
continuities,  would  slightly  shift  down  the  tunneling 
frequency,11  but  they  do  not  notably  affect  the  main  features 
in  this  scenario.  The  low  value  of  eeff  in  the  channel  implies 
/3ch  —  0 ,  which  ensures  small  phase  delay  and  uniform  field 
inside  the  channel,  whereas  the  continuity  of  voltage  across 
the  step  leads  to  an  enhancement  of  the  electric  field  inside 
the  channel  by  a  factor  of  1  /s  with  respect  to  the  field  in  the 
outer  waveguide  sections. 

Consider  now  a  small  cavity  region  with  permittivity 
ecav  A  ech  in  the  narrow  channel,  as  depicted  in  Fig.  1 .  Apply¬ 
ing  the  same  TL  model,  condition  (3)  is  modified  into 


y2ch\ 1  + 


2  sin[/3cavLcav](  yc,dW  Vch) 


i(^7ch  ^7cav  (  ^7ch  ^7cav)^^l-/^ch(X  ^ycav)])sifl[^cav^'cavl  "l"  ^  ^7cav^7ch  ^^^t^cav^cavl^i^C^chC^  ^cav)]} 


=  %  (4) 

s 
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DIELECTRIC  SENSING  IN  e-NEAR-ZERO... 

where  the  subscript  cav  refers  to  the  cavity  parameters.  For 
the  sake  of  simplicity,  but  without  loss  of  generality,  we  have 
assumed  that  the  cavity  is  positioned  in  the  center  of  the 
channel.  Equation  (4)  quantifies  the  perturbation  and  the  shift 
in  the  tunneling  frequency  caused  by  the  cavity.  In  particular, 
it  may  be  deduced  from  Eq.  (4)  that  the  shift  is  always  to¬ 
ward  lower  frequencies  if  ecav  >  ech,  implying  that  when  the 
cavity  is  inserted  the  tunneling  occurs  at  a  frequency  lower 
than  the  channel’s  cutoff  frequency.  It  is  interesting  to  point 
out  that  for  this  symmetric  scenario,  for  which  the  cavity  is 
placed  at  the  center  of  the  channel,  total  transmission  can 
still  be  achieved  when  Eq.  (4)  is  satisfied.  For  the  asymmet¬ 
ric  case,  however,  the  resonant  tunneling  does  not  ensure,  in 
general,  total  transmission. 

Figure  2  shows  the  transmission  coefficient  predicted  us¬ 
ing  the  analytical  TL  model  described  above  (dashed  line) 
and  compared  with  the  full-wave  numerical  simulations12 
(solid  lines)  for  different  geometries  of  the  channel  and 
the  cavity.  In  all  these  examples  L= 1 27  mm,  b=2a 
=  101.6  mm,  e=2e0  (the  permittivity  of  Teflon),  and  ech=e0- 
The  black  lines  refer  to  the  case  in  which  the  cavity  has  the 
same  permittivity  as  the  channel,  i.e.,  ecav=  ech=  e0,  which  is 
consistent  with  the  geometry  analyzed  in  Ref.  6.  The  ENZ- 
related  tunneling  occurs  very  close  to  the  cutoff  frequency  of 
the  channel,  which  happens  at /=  1.48  GHz  in  this  scenario. 
A  second  peak,  consistent  with  a  Fabry-Perot  resonance,  is 
obtained  at  a  higher  frequency,  but  this  is  strongly  dependent 
on  the  length  L. 

Figure  2(a)  considers  the  presence  of  a  cavity  with  Lcav 
=L/10  in  a  channel  with  ach=al 64.  Small  variations  in  ecav 
introduce  a  significant  shift  in  the  ENZ-related  tunneling  fre¬ 
quency,  in  accordance  with  Eq.  (4).  The  full-wave  numerical 
simulations  agree  very  well  with  the  analytical  model  de¬ 
rived  above  (which  neglects  the  presence  of  reactive  loads 
describing  the  steps  around  this  tunneling  frequency),  vali¬ 
dating  the  assumption  that  the  presence  of  an  effective  ENZ 
region  around  the  channel  E- plane  discontinuities  signifi¬ 
cantly  reduces  the  excitation  of  higher-order  evanescent 
modes.11  At  a  higher  frequency,  however,  when  the  Fabry- 
Perot  resonance  arises  (for  this  geometry  around  1.8  GHz) 
and  the  channel  has  a  larger  positive  eeff,  the  presence  of 
these  loads  cannot  be  safely  neglected,  as  revealed  in  the 
disagreement  between  the  analytical  and  simulation  lines 
around  this  frequency.  We  note  here  that  the  high-g  ENZ- 
related  transmission  peak  in  the  channel  may  be  fine  tuned 
by  varying  the  cavity  permittivity  ecav  around  the  tunneling 
frequency.  However,  the  Fabry-Perot  resonance  peak,  in 
comparison,  is  weakly  affected  by  the  presence  of  the  cavity 
due  to  the  nonuniform  distribution  of  electric  field  in  the 
channel  in  the  Fabry-Perot  scenario.  In  particular,  in  this  case 
the  electric-field  distribution  is  such  that  the  field  is  expected 
to  be  near  zero  in  the  center  of  the  channel  where  the  cavity 
is  positioned. 

Figure  2(b)  considers  a  similar  geometry  but  with  a  chan¬ 
nel  with  higher  height,  i.e.,  ach=al  16.  In  this  case  the 
electric-field  enhancement  in  the  channel  is  less,  and  corre¬ 
spondingly  the  Q  of  the  tunneling  transmission  is  also  low¬ 
ered,  causing  broadening  of  bandwidth  of  transmission 
peaks.  Figure  2(c)  presents  the  results  for  a  larger  cavity,  i.e.. 


FIG.  2.  (Color  online)  Transmission  coefficient  for  the  U  chan¬ 
nel  of  Fig.  1  with  L=127  mm,  b  =  2a=  101.6  mm,  e=2e0,  and 
ech=e0:  (a)  ach=a/64,  Lcav=L/10;  (b)  ach=a/16,  Lcm=LI  10;  and 
(c)  ach=fl/64,  Lcav=L/5.  Solid  lines:  full-wave  simulations  (Ref. 
12);  dashed  lines:  TL  model. 

Lch=LI 5.  In  this  case  the  resonance  shift  is  more  pro¬ 
nounced,  as  expected  from  the  larger  geometrical  perturba¬ 
tion,  but  maintaining  a  similar  Q  factor  as  in  the  case  of  Fig. 
2(a),  since  ach  remains  unchanged. 

These  results  are  summarized  in  Fig.  3,  which  reports  the 
calculated  shift  of  tunneling  frequency  [Fig.  3(a)]  and  the 
variation  of  the  transmission  coefficient  at  the  channel  cutoff 
frequency  [Fig.  3(b)]  for  the  three  geometries  of  Fig.  2.  This 


045102-3 


ANDREA  ALU  AND  NADER  ENGHETA 


PHYSICAL  REVIEW  B  78,  045102  (2008) 


FIG.  3.  (Color  online)  Variation  of  (a)  the  ENZ-related  tunneling 
frequency  and  (b)  the  ENZ-related  transmission  coefficient  at  the 
channel’s  cutoff  frequency  /=  1.48  GHz  vs  normalized  ecav  for  the 
geometries  of  Fig.  2. 

clearly  reveals  the  filtering  response  of  this  ENZ-related  tun¬ 
neling  phenomenon  and  suggests  that  the  ENZ-inspired  su¬ 
percoupling  may  be  exploited  as  a  sensor  for  permittivity 
measurement  by  properly  tailoring  the  tunneling  Q  factor. 
Operating  at  the  channel’s  cutoff  frequency,  the  transmission 
coefficient  can  be  sensitively  altered  when  a  small  object 
(particle,  fluids,  etc.)  with  permittivity  ecav  is  inserted  in  the 
channel.  By  measuring  this  transmission  coefficient,  one  can 
then  evaluate  the  value  of  ecav,  based  on  a  curve  such  as  Fig. 
3(b).  We  have  assumed  here  that  ecav  is  a  real  quantity,  since 
losses  in  dielectric  materials  may  be  safely  neglected  in  this 
context,  but  this  analysis  remains  valid  also  for  absorbing 
materials. 

Figure  4  shows  the  electric-  and  magnetic-field  distribu¬ 
tions,  normalized  to  the  impinging  field  amplitude,  calcu¬ 
lated  on  the  bottom  plate  of  the  waveguide  across  the  chan¬ 
nel  for  the  geometries  in  Figs.  2(a)  and  2(c)  with  ecav = 3 e0  at 
the  corresponding  resonant  frequencies,  as  indicated  in  the 
caption.  It  is  clear  that  we  have  a  large  increase  in  the  nor¬ 
malized  electric-field  distribution  inside  the  cavity,  even 
more  than  what  is  attainable  within  an  empty  supercoupling 
channel  at  its  cutoff  frequency,  which  for  these  geometries 
would  correspond  to  a  uniform  distribution  all  across  the 


FIG.  4.  (Color  online)  Normalized  (a)  electric-field  Ey  and  (b) 
magnetic-field  Hx  distributions  [real  (solid)  and  imaginary  (dashed) 
part]  on  the  bottom  plate  of  the  waveguide  at  the  tunneling  frequen¬ 
cies /=  1.327  GHz  and /=  1.206  GHz,  respectively,  for  the  geom¬ 
etries  of  Figs.  2(a)  and  2(c)  with  ecav=3£0.  The  channel  region  is 
highlighted  in  the  plots. 

channel  with  normalized  amplitude  equal  to  ,v  =  64.  With  the 
cavity  included,  however,  the  electric  field  is  further  in¬ 
creased  in  the  central  region  of  the  channel,  where  the  cavity 
is  placed,  due  to  an  anomalous  growing  exponential  tail  in 
the  first  part  of  the  channel  associated  with  resonant  tunnel¬ 
ing  and  the  fact  that  the  channel  is  resonating  below  its  cut¬ 
off.  This  is  even  more  evident  in  the  magnetic-field  distribu¬ 
tion  [Fig.  4(b)],  which  would  be  constant  in  an  empty  ENZ 
channel.4  The  downshift  of  the  tunneling  frequency  produces 
a  large  buildup  of  the  normalized  electromagnetic  fields  in 
the  channel,  localized  at  the  cavity  edges.  In  some  sense,  this 
is  analogous  to  the  exponential  growth  of  the  electromag¬ 
netic  fields  in  the  resonant  pairs  of  metamaterials  below 
cutoff,13  even  if  in  this  scenario  the  resonance  interestingly 
happens  in  a  simple  hollow  waveguide  section.  The  field 
increase  is  more  pronounced  for  a  larger  cavity,  since  the 
resonance  happens  at  lower  frequencies  for  which  the  effec¬ 
tive  permittivity  of  the  channel  is  more  negative,  following 
Eq.  (1). 
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HI.  CONCLUSIONS 

We  have  reported  here  an  idea  to  exploit  the  dramatic 
field  enhancement  associated  by  energy  squeezing  and  tun¬ 
neling  in  metamaterial-inspired  ultranarrow  waveguide  chan¬ 
nels  with  near-zero  effective  permittivity  in  order  to  sense 
small  permittivity  variations  in  a  tiny  object.  We  have  ana¬ 
lytically  modeled  the  supercoupling  effect,  deriving  closed- 
form  expressions  to  describe  the  perturbation  caused  by  pres¬ 
ence  of  defects  or  permittivity  perturbations  along  the 
channel.  The  experimental  microwave  realization  of  dielec¬ 


tric  sensing  presented  in  this  paper  is  currently  under  test  in 
our  group,  and  the  extension  of  these  concepts  to  optical 
frequencies  is  also  conceptually  viable  using  plasmonic 
waveguide  channels.  Applications  for  tuning  the  supercou¬ 
pling  pass-band  frequency  and  for  h  i g h - Q  sensing  are  envi¬ 
sioned  in  different  frequency  regimes 
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